A Complex Chemical Potential: Signature of Decay in a Bose-Einstein Condensate 
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We explore the zero-temperature statics of an atomic Bose-Einstein condensate in which a Fes- 
hbach resonance creates a coupling to a second condensate component of quasi-bound molecules. 
Using a variational procedure to find the equation of state, the appearance of this binding is mani- 
fest in a collapsing ground state, where only the molecular condensate is present up to some critical 
density. Further, an excited state is seen to reproduce the usual low-density atomic condensate 
behavior in this system, but the molecular component is found to produce an underlying decay, 
quantified by the imaginary part of the chemical potential. Most importantly, the unique decay rate 
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) can be measured in experimental 



Since the production of atomic Bose-Einstein conden- 
sates (BEG) in the laboratory niany schemes have 
been proposed whereby the experimentalist may control 
the interatomic interactions governing the behavior of 
these gases One such proposal involves a Feshbach 
resonance in which two atoms combine to form a quasi- 
bound molecule 0, |^ . This molecule is described as the 
intermediate state or closed channel of the scattering re- 
action as the constituent atoms in general have differ- 
ent spin configurations in the bound state than in the 
scattering state. Due to its dependence on the internal 
spin states, the energy difference, or detuning between 
the scattered and bound states can thus be tuned using 
the Zeeman effect in an external magnetic field. As the 
binding energy of the molecular state is brought close 
to the energy of the colliding atoms, the appearance of 
these loosely bound molecules increases. Gonsequently, 
the coupling between atoms and molecules acts to modify 
the effective interatomic interactions. Near zero energy, 
these interactions are described by the s-wave scattering 
length, which in turn can be tuned by varying the exter- 
nal magnetic field. This degree of control suggests that 
the negative scattering length of initially unstable con- 
densates may be tuned to positive values thereby render- 
ing the condensate stable. We provide a many-body vari- 
ational description of such a scenario, specifically using 
the case of ^^Rb to compare our results with experiment 
0, Q • For a uniform system our findings reveal a collaps- 
ing ground state and a decaying excited state, where the 
latter fits the behavior seen from experiment. In the ex- 
cited state, a complex chemical potential is obtained in 
which the imaginary part determines the inverse decay 
time. 

Before embarking on the many-body analysis of the 
coupled atom- molecule BEG, it is first necessary to re- 
view the relevant two-body physics underlying the inter- 
atomic interactions. It can be seen that the two-body for- 
malism provides a means by which the many-body equa- 
tions can be renormalized by replacing the interaction 



strength by the s-wave scattering length as the relevant 
parameter. In two-atom scattering, there are a number 
of different channels or outcomes of the scattering event, 
corresponding to both the closed and open channels of 
the bound and scattering states, respectively. Because 
coupling will be strongest to the highest lying molecular 
state below the continuum of scattering states, we con- 
sider a two-channel model with coupling between a single 
open and a single closed channel. Since the coupled chan- 
nels analysis has been described in detail elsewhere 0], 
we present only the main results of our approach. Using 
a separable form for the the two-body pseudopotential 
we have, 



(k|F|k') = Aw(k)t;(k') 



(1) 



where k is the momentum in the center of mass and A 
is the interaction strength which is taken to be negative, 
since we assume an attractive background throughout j3| . 
Using the separable potential with a molecular form fac- 
tor equal to w(k), the coupled channels analysis yields an 
integral equation for the two-body wavefunction, ^'(k), 



(2k2 - *(k)+ {\ - j j <?^fc'w(k')^(k') = 0, 

(2) 

where E, A and e are the center of mass energy, in- 
teraction strength and detuning, respectively. Here 
and throughout all energies are on the scale of /2m. 
Equation Q is recognized as the usual single chan- 
nel two-body integral scattering equation where the 
molecular coupling has replaced the background inter- 
action strength. A, with an energy dependent strength, 
X — /^a^ / {e — E). One striking feature of this equation 
lies in the observation that for some nonzero coupling, 
g, the interatomic interaction may be produced by a sin- 
gle Feshbach resonance, for which the scattering equation 
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becomes 



(2k2 - E)-^{k) = v{k) J d^k'v{k')-^{k'). (3) 

If the detuning, e, is far enough away from the two-body 
energy, | e| ^ | then the effective interaction strength 
is independent of energy {—g^/ie — E) « —g^/e) and can 
be chosen to be equal to A, the strength of the original 
pseudopotential. Thus, the two-body background poten- 
tial can be modeled by a single Feshbach resonance with a 
detuning that is very far from the system's energy. It will 
be seen that this form of the interatomic force results in 
a simplified expression for the two-body interaction op- 
erator in the many-body Hamiltonian. 

Starting with the relationship between the scattering 
length, a, and the T matrix, Stto = (k|T|k') |k=k'=o, we 
can derive the renormalization equation for the interac- 
tion strength A. The separable potential of Eq.Q indi- 
cates a separable form of the T matrix as well, (k|T|k') = 
vik)tv{k'). Thus, using the Lippmann-Schwinger equa- 
tion, T = V + VGT, where G is the free particle Green's 
function, we can solve for the T matrix and, consequently, 
the background scattering length, abg. Subsequently, we 
take the zero range limit for the interatomic potential in 
which we allow v{k) — > 1, thus giving 



a ^ Q/37 

+ -^Y.^xiF^Pyi^^i^p + H.c.). (6) 
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In anticipation of the variational procedure to be car- 
ried out, a term ~fiN has been added where the Lagrange 
multiplier is identified as the chemical potential. Further- 
more, the atomic field is coupled to both the molecular 
field, (j), and the auxiliary field, x, with coupling con- 
stants {Xa,g) and detunings (e,£) defined by Eq.© and 
Eq.lPJ, respectively. Since the zero range limit is taken, 
we can, without loss of generality, take the same molec- 
ular form factor, F, setting it equal to the v defined in 
the separable potential, Eq.lQJ. As discussed, the back- 
ground atom-atom interaction is reproduced by setting 
5^/(£ — E) = —A with \ e\ :s> \ E\. In the zero range 
limit (A — > 0^), these conditions are satisfied if we set 
e = 2/A2 and g = ^/-2/X. 

For the zero-tcmpcrature analysis of the Hamiltonian 
|l3l |. we work in the Schrodinger picture using a static, 
" squeezed" Gaussian trial wave functional given by 
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(4) 



where the range, 6, tends to zero as A — > 0^ in such 
a manner that atg remains negative and finite |inl |. As 
found from Eq.(|2Jl, we replace the background strength 
A by A — }^o? I (e — E) to include the effect of molecular 
coupling. With _B = 0, we use this substitution to obtain 
the full scattering length 



aie) 



eabg 



(5) 



Alternatively, the full scattering length can be expressed 
in terms of the external magnetic field, B, since the de- 
tuning is given as e = ^{B — Bl), with 7 being the atomic 
species dependent proportionality and with Bl as the off- 
set field. These two-body results will be used to renor- 
malize the many-body treatment. 

In the many-body analysis, we begin with a Hamilto- 
nian which includes the atomic kinetic energy, atom-atom 
collisions, the molecular detuning energy and a coupling 
between atoms and molecules. However, according to 
the discussion behind Eq.Q, we can replace the usual 
two-body term, t/j'^, with an auxiliary field, coupled 
to the atomic field, ip, through a single Feshbach reso- 
nance, thereby giving a simplified Hamiltonian which is 
quadratic in the atomic field instead of quartic [ill : 
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N^exp I - E } ^xe^P | " E ^-^^^^ } ' (^) 

where N^, and N^. are normalization constants, and 
the width, G, is taken as real. Denoting -0' as the atomic 
field with "0 as its mean, the atomic field fluctuations are 
given as (5-0' ^ tp' — ipj with the analogs of the molecular 
fluctuations defined similarly. In the momentum space 
representation of the uniform system, all fields assume 
their mean values times a Dirac delta function, S{k), 
whereas the Gaussian width assumes the diagonal form, 
G(k)~^. Taking the expectation value of Eq. © with 
respect to the trial functional of Eq. (TJ gives an expres- 
sion for the pressure, —P = {H — fiN) /V, in a volume 
V. Extremizing the pressure in G(k), and in the mean 
fields, ip, (f>, and x, gives a set of variational equations, 
which in momentum space can be expressed as 



1 k^-^, + p 



(/3-m)^ = 0, 



(8) 



(9) 
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d^kD{k) + -v^ 



0, (10) 



where we have defined f3 = a(f> + gx- Additionally, we 
denote the atomic density, {■ip^ip)/V, and the anomalous 
atomic density, {■tp'4!)/V, by /i?(k) +-02/2 and jD{k) + 
0^/2, respectively, where the fluctuations -R(k) and D(k) 
are related through the Gaussian width as 



i?(k) = ^ 



iG-i(k) + G(k)-l 



^(k) - ^ 



-G-i(k)-G(k) 



(11) 



(12) 



Using Eqs. (|S|)- ((T^ . expressions for the energy density, 
u — (H) /V, and the mass density, p, are obtained. 



■0^/2. In this case, the energy per particle starts at pc, 
then continues downward for arbitrarily large p, (Fig. 

Thus, these two solutions comprise a single curve of 
two pieces, recognized as the collapsing ground state of 
the model Hamiltonian. 



> 

CD 



CO 
Q. 



0.5 



-0.5 



\|/7t0 




Density (cm ) 



xlO' 



-f^-^+/?y<^^fcw, (13) 

p = /.3,,,(k) + lv.^ + ^^. (14) 

In the zero range (6 — > 0,A — > 0^) limit, it can be 
seen from Eqs. (|ll|l and 112() that the divergent terms in 
the integrals in u have the behavior Jd^kD(k) ~ (3 /'2b 
and Jd^kk'^Rik) ~ (3^ /4b. Using Eq. Q relating A 
to the range, b, the divergences cancel, resulting in an 
expression dependent only on abg. This renormalization 
procedure is used in all solutions given by Eqs. 

We begin the analysis of these solutions by noticing 
that for a nonzero P in Eq. (jSJ, /i must be negative in 
order to keep G real. Taking a negative p, we see from 
Eq. lO that there are two solutions, with one given by 
a vanishing atomic field (0 = 0), and the other given by 
a non- vanishing field and (3 — p < 0. Considering the 
■0 = solution first, Eqs. Hll|) and Q are substituted 
into Eq. (|10() , resulting in an expression that determines 
p in terms of (3. We find that this equation has a 
solution for f3 from zero only up to some critical value, 
Pc- From Eq. 1)14(1 . this range of beta corresponds to a 
range of densities, starting from zero with P and ending 
at some critical value, pc, corresponding to Pc- As shown 
in Fig. ^ the plot of the energy per particle, e — u/p, 
for this solution results in a downward curve starting at 
half the molecular binding energy, BE/2, and stopping 
at Pc- To find the remaining solution for p > pc, we take 
P = p in Eq. ©, and then use Eq. ^fTU\i to solve for 



FIG. 1: The real part of the excited state (a) and the collaps- 
ing ground state (b) are evaluated using the following numer- 
ical values for *^Rb : applied magnetic field, B = 162.3 G 
offset field, Bl — 165.6 G, coupling constant, —SnalgO^ = 
330, detuning proportionality, 7 = —30 G~^, and the respec- 
tive background and full scattering lengths, abg = —450 ao 
and a(e) ~ 193 ao, where ao is the Bohr radius p3 . With 
these values, we obtain a critical density, pc, of 1.22 x 10^^ 
cm~'^ and a binding energy, BE, of —3.7 neV. As the applied 
magnetic field is decreased toward resonance, pc decreases, 
thereby shortening the range of the = solution. Tuning 
the scattering length to negative values results in a collaps- 
ing atomic-molecular condensate solution, with the usual low 
density dependence, e ~ 47ra(e)p. 

Thus far, we have not uncovered the solution cor- 
responding to the low density behavior, e ~ 47ra(e)p, 
reported in the experiments for the positive scattering 
length regime of ^^Rb ,6]. Obtaining this solution re- 
quires a closer examination of the density expression for 
P ^ p. Using Eqs. ©, l|inil-l|I3), the right hand side 
of Eq. 114|) can be expressed in terms of the chemical 
potential, p, 

(15) 

Since Eq. (|15|l has no solution for both p, p > 0, we 
expand in p, which is constrained to be real and positive, 

p = 8TTa{e)p - i^256a(e)^/2p3/2 + . . . ^ (16) 

Thus, this expansion is seen to give the desired low- 
density 47ra(e)p dependence of the energy per particle, 
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but becomes complex at higher order [l5j. 

At first glance, the complex-valued chemical potential 
appears unphysical. However, the chemical potential can 
also play the role of phase of the atomic field, "0) as can 
most easily be seen from the Heisenberg equation for ip, 
ihdil^*/dt = {[^\H]) = -I3{t)*ip{t), in which the expec- 
tation value is carried out using the trial wave-functional 
of Eq. IJj). Taking the time dependence of the fields as 
tp{t) = tpexp{—ifit/h) and P{t) = Pexp {—2ifit/h), yields 
the variational Tp equation Q. 

From this vantage point, the imaginary part of fj, is seen 
to be proportional to the inverse of the decay constant. 
The inspiration for such an interpretation arises in the 
context of QED, where a complex action signifies decay 
of a constant, uniform electric field Using the ^^Rb 
parameters given in Fig. ^ and taking ^ 10* atoms in 
a BEC cloud of radius 25 /im 6], a decay constant, t, 
of approximately 14 s is obtained to lowest order in p 
(r ^ 1/p^/^). The fuU solution of Equation (dJ leads to 
a seventh order equation in y/~iJ-, which, when solved nu- 
merically, gives a decay constant of approximately 14.2 
s, thus confirming the validity of the expansion in Eq. 
(fTB|l . It must be pointed out that this is a novel decay 
process inherent in the Hamiltonian describing only the 
two-body interactions in the gas. This result is in quali- 
tative agreement with the 10 s decay reported in Ref. 0, 
where the experiments were carried out in a regime where 
this decay process dominates over the two and three-body 
inelastic processes [l^l • As seen from Eq. ((TB)l , this decay 
rate has novel density and scattering length dependencies 
a(e)^/^p'^/^) which can be experimentally tested. 

In summary, we have used a variational procedure 
in extremizing a many-body Hamiltonian describing an 
atom-molecule BEC coupled through a Feshbach reso- 
nance. For positive scattering lengths, we have found 
that the existence of the bound molecular state results 
in a collapsing ground state and an excited state with 
a complex- valued chemical potential. By examining the 
Heisenberg equation for the atomic field, the imaginary 
part of the chemical potential is seen to physically cor- 
respond to the inverse decay time. Finally, this decay 
process was found to give the dominant contribution in 
sufficiently dilute systems. 
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